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Abstract 

Within the continuum dislocation theory the asymptotic analysis of the plane 
strain crack problem for a single crystal having only one active slip system 
on each half-plane is provided. The results of this asymptotic analysis show 
that the square root stress singularity remains valid during the plastic defor¬ 
mation, while the dislocation density is proportional to the stress intensity 
factor and distributed as the square root of the distance from the crack tip. 
The analytical solution for the angular distribution of the dislocation density 
is found. 

Keywords: crack, dislocations, plastic slip, stress singularity, dislocation 
density. 


Dislocations appear to reduce energy of crystals. For crystals with cracks 
the high stress concentration near the crack tip causes also high energy of 
crystals in that region. It is therefore natural to expect that, when the load 
is sufficiently large, dislocations nucleate near the crack tip to reduce the 
stress level and by this also the energy of crystals. It is then crucial to have 
the correct perception of how dislocations nucleate near the crack tip. Up 
to now, the commonly accepted point of view is that dislocations nucleate 
directly at the crack tip and then glide away from it under the Peach-Koehler 
force m- However, the analysis of crack problems reveals that the resolved 
shear stress is large not only at the crack tip, but also in its neighborhood. 
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Taken this for granted, then, according to the Schmid’s law, dislocations 
must appear simultaneously in that neighborhood exhibiting the collective 
character of dislocation nucleation. Since the typical dislocation density is 
high (about 10 8 -t-10 15 dislocations per square meter), it makes sense to use 
the continuum approach to study this problem. 

This short paper aims at finding the stress and dislocation distribution 
near the crack tip in ductile single crystals within the continuum dislocation 
theory (CDT) proposed by Berdichevsky [3 0] and developed further in [5]- 
pL5 j . Considering the plane strain crack problem, we assume that during the 
plastic deformation only one slip system on each half-plane of the crystal is 
active. We provide an asymptotic analysis of this crack problem in the polar 
coordinates. The results of this asymptotic analysis show that the square 
root singularity for the stress field near the crack tip remains valid. This 
agrees with the singularity of HRR-keld obtained by Hutchinson [16], Rice 
and Rosengreen in conventional plasticity for the materials with linear 
hardening. What the dislocation distribution near the crack tip is concerned, 
we show that they must be distributed such that the resolved shear stress is 
balanced with the back stress in accordance with the equilibrium of micro¬ 
forces acting on dislocations. This leads to the power law distribution y/r, 
with r being the distance from the crack tip, for the dislocation density, 
with its intensity being proportional to the stress intensity factor. We find 
also the universal angular distribution of the dislocation density. Note that 
the crack-tip fields in single crystal has been analyzed within the discrete 
dislocation dynamics in |18j . Experimental observations of the dislocation 
distribution near the crack tip in single crystals by electron tomography 
have been reported in [T5]- Another quite promising experimental method 
of measuring the dislocation density by using electron backscatter diffraction 
(EBSD) technique has been developed in [20.. 2T] . 

1. Plane strain crack problem for single crystal 

Consider the plane strain problem for a single crystal containing a crack 
lying on the left-half of the aq-axis as shown in Fig. [TJ The depth of the crystal 
in the aq-direction is taken large enough to guarantee the plane strain state 
having two non-zero components of displacement vector U\ = and 

«2 = U 2 (xi,X 2 ). The crystal is oriented in such a way that its lattice and 
mechanical properties as well as the loading condition (say tractions acting 
at the outer boundary) are symmetric with respect to the reflection about 
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Figure 1: The plane strain crack problem for single crystal 

the ay-axis. Because of this mirror symmetry it is sufficient to consider the 
upper-half of the crystal. If the load is small enough, then it is natural to 
assume that the crystal with this crack deforms elastically. However, if the 
load exceeds some threshold value, dislocations can occur causing the plastic 
deformation of the crystal. We assume that, during this plastic deformation, 
only one slip system from each half of the crystal is active and the dislocations 
are straight lines parallel to the ^ 3 -axis. For the more realistic crack problems 
in single fee and bcc crystals having several active slip systems the reader 
may consult [ 22 ]- Letting s = (cos ip, sin (p, 0) denote the slip directions, 
and m = (— sin ip, cos (p, 0) the normal vector to the slip planes of the slip 
system in the upper-half of the crystal, we may express the plain strain plastic 
distortion tensor in the form /3 = /3{x\, X 2 )s <8) m. We are going to determine 
the displacements Ui(xi,X 2 ), u^Xi, ^ 2 ), and the plastic slip f3(xi,x 2 ) near 
the crack tip during this plastic deformation. 

For the plane strain state the non-zero in-plane components of the sym¬ 
metric strain tensor e = |(Vu + uV) are 

£11 = W,i, £12 = £21 = 2^ Ul ’ 2 M 2 >i)> £22 = u 2 ’ 2- 

Throughout the paper the comma standing before an index is used to denote 
the partial derivative with respect to the corresponding coordinate. The in¬ 
plane components of the symmetric plastic strain tensor e v = |(/3 + j3 T ) 
equal 

£11 = ~\P sin 2< /b £ i2 = £ 2i = cos 2( P, e 22 = sin 2i P- 
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With these total and plastic strain tensors we obtain the in-plane components 
of the symmetric elastic strain tensor s e = e — e p in the form 


£n = u 1,1 + ~/3 sin 2ip, £ 22 = ^ 2,2 - x/3 sin 2 ip 

2 ! ' (1) 
£ 12 = £ 21 = 2^ Ul ’ 2 + U2J “ ^ C ° S 

Let us compute Nye-Bilby-Kroner’s dislocation density tensor (introduced in 
[23j-[25j) ol = —/3 x V, with x being the vector product. For the plane strain 
plastic slip f3(xi,x 2 ) there are two non-zero components of this tensor given 
by 


op 3 = (Ai cos ¥ + /3,2 sin ip) cos ip, a 2 3 = (/?, 1 cos ip + /3 2 sin p>) sin ip. 

These are the component of the net Burgers’ vector of all excess dislocations 
whose dislocation lines cut the area perpendicular to the x 3 -axis. Thus, 
the net Burgers’ vector of excess dislocations shows in the slip direction s 
indicating that we are dealing with the edge dislocations only. The scalar 
dislocation density (or the number of dislocations per unit area) equals 

P= \\/ a i3 + "23 = r\P,iCOS(p + /3 2 sin</?| = J|d s /3|, (2) 

where d s (3 = V/3-s. Similar quantities in the lower-half plane can be obtained 
if we replace in the above formulas ip by — p. 

For most metals the elastic strain tensor s e is usually small. Therefore 
the free energy density per unit volume of the crystal with continuously 
distributed dislocations can be proposed in the form M 

p) = 7^(>£ e ) 2 + P £e -£ e + l lk In 

Ps 

with A and /j the Lame’s constants, p s the saturated dislocation density, and 
k the material constant. The first two terms in (|3| represents the contribution 
to the energy due to the elastic strain, while the last term corresponds to the 
energy of the dislocation network. The logarithmic energy term stems from 
two facts: (i) energy of the dislocation network for small dislocation densities 
is the sum of energy of non-interacting dislocations (see the reasonings based 
on the statistical mechanics of dislocations in [5B; [2T], and (ii) there exists 
a saturated dislocation density which characterizes the closest packing of 
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dislocations admissible in the discrete crystal lattice. The logarithmic term 
[4] ensures a linear increase of the energy for small dislocation density p 
and tends to infinity as p approaches the saturated dislocation density p s 
hence providing an energetic barrier against over-saturation. For small up to 
moderate dislocation densities this logarithmic term can be replaced by its 
Taylor expansion 
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For simplicity of the subsequent analysis we shall use this approximation only. 
Taking all the above formulas into account, we write the bulk energy density 
per unit volume of the crystal with continuously distributed dislocations in 
the form 


-0(s e , p) = ^A(wi,i + u 2)2 ) 2 + p{ui ,i + sin 2ip) 2 + p(u 2 p ~ sin 2<p) 2 

,i, , * 0 , 2 . , 1 (d s p)\ 

+ -P( u l,2 + U 2,l — P cos 2(p) +/ifc(-^—-f - ). 

The energy functional per unit depth of the crystal becomes 


I[u a ,f3] = 


7 } A(wi.i + U 2 , 2 ) 2 + p(ui,i + sin 2 cp) 2 + p,(u 2;2 - ^ f3 sin 2ip) 2 


1 


+ 1,2 + Uo i - f3 COS 2<p) 2 + /i/c( 

2 op, 


W\ , i(a s /3) 2 


+ X 


2 Vpl 


dx, (4) 


where hi is the domain occupied by the undeformed crystal having the crack 
and dx = dx\dx 2 . Since we are interested in the fields near the crack tip, 
the virtual work of the traction acting at the outer boundary is dropped 
in this functional. We require that the displacements U\ and u 2 and the 
plastic slip /3 be continuous on the positive aq-axis. On the contrary, on the 
crack faces x\ < 0, x 2 = ±0 (which are assumed to be not in contact with 
each other in the deformed state for the mode I crack) no constraints are 
imposed on U\, u 2 , and /3 (free boundary). Provided the dissipation caused 
by the dislocation motion is negligible, then the true displacements u±, u 2 
and plastic slip j3 in the final state of equilibrium minimize energy functional 
Q among all admissible displacements and plastic slips satisfying the above 
boundary conditions. 

Let us derive the equilibrium equations and boundary conditions for the 
true displacement vector and the plastic slip. For plane strain problems it is 
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convenient to use the index notation with the greek indices running from 1 
to 2 and with the summation convention being applied for two-dimensional 
components of vectors and tensors. Then, for instance, the elastic strain 
tensor has the following non-zero components 

4/3 = + Ufa) - ^ 3(s a mp + spm a ). 

Keeping this index notation in mind, we compute now the first variation of 
the energy functional (|4]) 


51 = 


o-, 


*p( 6 u aj i 3 - 5/3s a mp) + K a 5/3^\ dx, 


where 


<9-0 e 

°a/3 — a c e — + ZfJ,£ a p, 


Krv = 


90 


= pk 


de U 5 

1 ' (A R\ . ds @ 

—sign (d s p) + — 

bps b-p 2 


(5) 


We call <j a p the (symmetric) Cauchy stress tensor, and n a the higher order 
stress vector. Integrating by parts and omitting the surface integral at the 
outer boundary, we obtain 


61 



(r + K a ,a)6/3] dx 


' xi<0 


(a a 2 Su a + K 2 6 P)\l 2 2 Jlodxi, 


where r = cr a ps a 7 n /3 is the resolved shear stress (or Schmid stress). Thus, the 
equation 51 = 0, together with the arbitrariness of 5u a and 5/3, implies the 
equilibrium equations of macro-forces acting on the material volume element 


^"a/3,/3 0) (6) 

and of micro-forces acting on dislocations 


T + K a ,a — 0 


(7) 


in as well as 


cr a2 = 0, k 2 = 0 for x\ < 0, x 2 = ±0. 


( 8 ) 


6 







In addition to these we can pose the boundary condition 


f3(xi, 0) = 0 for x\ > 0 


( 9 ) 


which is simply the consequence of the mirror symmetry of the crack problem. 
The system of equations ([5]), ([6]), ([7| together with the boundary conditions 
(|8]) and (J9]) constitute the local crack problem to determine the displacements 
u a and the plastic slip /?. 

To get the governing equations in terms of displacements and plastic slip 
we substitute the constitutive equations ( |5|), with e e a p and p from (JTJ) and 
(|2]), into the equilibrium equations (J6]) and ([7]). This leads to 

(A + p)up^a + iiUafip - n/3,p(s a mp + spm a ) = 0, (10) 


and 


T U/3 t a^)S a 7Tlp pf5 T — 2 ,ot/3^a^/3 O’ 


b 2 P. 


( 11 ) 


2. Asymptotic analysis 

Near the crack tip different terms in the governing equations (10) and ( fllj ) 
will have different orders of smallness. It is convenient to do the asymptotic 
analysis of these equations in the polar coordinate system r and 6 , where 


X\ — r cos 9, x 2 = r sin 9. 


Let us assume the asymptotically main terms for displacements and plastic 
slip near the crack tip in the form 


u a (r,9) = r n f a (9), (3(r,9) = r m g(9), (12) 


where f a (9) and g(9) are unknown functions describing the angular distribu¬ 
tion of the displacements and plastic slip, and n and m are unknown numbers. 
Since the displacements and the plastic slip cannot be singular at the crack 
tip, both m and n must be positive. Besides, based on the solutions of other 
similar crack problems we may assume that n lies in the interval (0,1). We 
substitute asymptotic formulas (12) into the governing equations and use the 
transformation rules 


3i ' 


cos 9 

— - sin 9 
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d r 

d 2 


sin 9 

- cos 9 
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. d ° . 
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to compute the derivatives. According to these transformation rules, the first 
derivatives reduce the power of r by one, and the second derivatives by two. 
Thus, the asymptotically main terms in (10) are the first two because the 
last one has the order r m_1 which is small compared with r n_2 . Neglecting 
this last term in equations (10) we get Lame’s equations of linear elasticity 
that are uncoupled from equation (11). Likewise, the boundary conditions 
after neglecting the last small term containing /3, reduces to 

\Ury^ 6 a 2 + H(u a! 2 + «2,a) = 0, 


so we get the plane strain elastic crack problem that can be solved by using 
the Airy stress function. The well-known solution of this problem shows that 
n — 1/2 (see, e.g., [27]). Besides, for mode I crack 


Kj nr e ( 2 e\ 

/i V 2tt 2 V 2 ) 

AT/ rr e ( 2 e\ 

H V 2tt 2 V 2 J 


1 


and 


cr ii 

0^22 


Kf 

y/2irr 

Kj 

y/2irr 


cos 2 


cos 2 


9 3 9 

1 — Sin - Sill — 
2 2 _ 

9 39' 

1 + sm - sm — 
2 2 


Kj 9 9 39 

a 12 = . cos - sm - cos —. 

pm 2 2 2 


(13) 


Here Kj is the stress intensity factor that can be found only after solving 
the global crack problem. Since crystals with continuously distributed dis¬ 
locations exhibit a linear work hardening as shown in [5]-|10j. this square 
root stress singularity agrees with the classical result obtained by Hutchin¬ 
son [16], Rice and Rosengreen H3 in conventional plasticity. Note also that 
the account of dissipation as proposed in [5]- [TO] does not affect this result. 

We turn now to equation (|TTj) . If m 6 (0,1], then the asymptotically 
principal term in this equation is obviously the last one, because it has the 
order r m-2 which is much larger than the orders of the first two. Neglecting 
the small terms in ( [TT| ) we get a simple equation 

b 2 p 2 s ds^ 
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The continuous solution of this equation 

(3 = bp 0 s , 


with s being the distance from the aq-axis in the s-direction and p 0 a constant 
dislocation density, would satisfy the equation and the boundary condition 
([9]) identically, but violate the boundary condition 


k 2 = pk 


—sign(<9 s /3) + 
bps 


dsl 

b 2 P 2 s . 


s 2 = 0 


on the crack faces. This is not surprising, because we know that the free 
boundary attracts dislocations, so the constant dislocation distribution can¬ 
not stay in equilibrium near the crack faces. Therefore m cannot lie in the 
segment (0,1]. For m G (1,2) the back-stress must be balanced with the 
resolved shear stress (having the power r -1 / 2 ) to guarantee the equilibrium 
of micro-forces acting on dislocations. Therefore m = 3/2 and (3 = r 3 ^ 2 g(6). 
Now we substitute this Ansatz into equation ( [TT| and maintain the terms of 
order r " 1 / 2 in it. Then equation ( [II] ) becomes 

/ y n 

OapScPB/j + (/? llSl + 2(3 ,\ 2 S\S 2 + /3,22^2) = 0; (14) 


where the stress components must be taken from (13). To compute the back 
stress we use the following formulas 


An = r-V» \ l 


3 1 

-cos 2d 

2 2 


M ~ ^ sin 2 dg\6) + sin 2 (%"(d) \ , 


A12 = r - 1 ' 2 
( 3,22 = r- 1/2 


- sin 2 6g(6) + - cos 2 6g'(0) — - sin 2 6g"(6) 
8 2 2 


3 1 

—I— cos 26 
2 2 


9(0) + 2 sin 2 0g\6) + cos 2 {6)g'\6) 


Substituting the resolved shear stress and back stress into (14) and simplify¬ 
ing this equation we obtain 


b 2 p 2 K r 

- s -== cos(3d/2 — 2<p) sind + 

pky/2'K 


3 

2 


3 

2 


1 

2 


cos(2d — 2 <p) 


9(0) 


— sin(2d — 2 tp)g'(d) + 2sin 2 (d — < p)g"(9 ) 


0. 
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Since this equation is linear, its solution g{9) must be proportional to the 
factor b 2 p 2 Ki / (yk\/2n). Therefore the problem reduces to solving the equa¬ 
tion 


cos(30/2 — 2ip) sin 9 + - 


cos ( 29 - 2( d) 


h{9) 


sin(26 l — 2 ip)h'{6) + 2sin 2 (6 ) — <.p)h"{6) = 0. (15) 


The inhomogeneous equation (15) must be subjected to the boundary con¬ 
ditions 

3 

h(0) — 0, - coscp h(n) + sin iph!(ii) — 0. (16) 

The second boundary condition means nothing else but the vanishing dis¬ 
location density near the crack faces. Only in this case the equilibrium of 
dislocations is guaranteed. 

To solve equation ( [l5| ) we make the change of unknown function 
h{9) = \/sin(0 — tp) y(9) for 9 > <p. 


It is straightforward to show that equation (15), in terms of y{6 ), becomes 

// , cos (30/2 -2ip) sin 9 

y + y = ~ 2 (s in (6 *-pw ' ( ) 

The general solution of ( Jl7| ) that can directly be verified, reads 
1 f) S/9 S/9 

y(9 ) =- = [—3sin - + sin ——4sin( ——2<^)]+Ci cos 9 + C 2 sin 9. 

3^/sm(9 — ip) 2 2 2 


For 9 < ip the change of unknown function h{9) = — \J sin((^ — 9) y{9 ) would 
do the job in that interval. Then we must combine the solutions in two 
intervals taking into account the continuity of the dislocation density and 
the analytic continuation of f3 across the singular point 9 = ip. Returning to 
the original function h{9) we obtain 


h(9) 


p(9 ) + (Ci cos 9 + C 2 sin 9) -y/sin(0 — p>) for 9 > ip, 
p{9) — (Ci cos 9 + C 2 sin 9) \J sin(<^ — 9) for 9 < <p, 


where 


p(9) 


l r „ . 9 .39 A . .39 

-1-3 sm - + sm — - 4sm(y 


2^)]- 
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Figure 2: Angular distribution of the plastic slip 


The unknown coefficients C\ and C ’2 can be found from the boundary 

4 sin 2p 


conditions (16) yielding 
Ci = 


3\/sin^’ 
Thus, the final solution reads 


4 sin 2 ip 

C 2 = ~ r . - cot p. 

o Y SHI ip 


m = 


p(d) 


4 sin 


i=(cos0 — cot p sin0)-y/sin(0 — p) for 6 > p, 


3^/sini^~ r > V V” t / ” r ’ o\ 

4sm =E(cos0 — cot p sin 0)y^sin(</?■— 6) for 6 < p. 


p{9) + , r .— 

t ^- rv ' oysin ip 

The plot of function h{9) (for p = vr/6) is shown in Fig. [2j 



Figure 3: Angular distribution of the dislocation density 
To summarize, the asymptotic formula for j3(r, 9) is 
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with h{9) from (18). Differentiating this formula, we find the signed disloca¬ 
tion density in the simple form 


P = 9 -f = ((0 - 2*0/2). 


(19) 


We see that the dislocation density is proportional to the stress intensity 
factor and distributed as yfr in the radial direction. Fig. [ 3 ] (again for <p = 
7t/ 6) shows the angular distribution q{6) = 2 cos 2 (0/2) sin((0 — 2yj>)/2) of the 
dislocation density. 


3. Conclusion 


We have shown in this paper that, within the continuum dislocation the¬ 
ory, the crack in ductile single crystals causes the square root stress singu¬ 
larity even during the plastic deformation. The dislocation density is pro¬ 
portional to the stress intensity factor and distributed in accordance with 


formula (19). The near-crack-tip fields for the displacements and plastic slip 


can be used to design singular crack-tip elements to solve crack problems 
in ductile single crystals by the finite element method. It remains still un¬ 
clear, whether some threshold in loading exists for the onset of dislocations 
nucleation. To answer this question one needs to compare the energy of the 
crystal containing a crack without and with dislocations. Another open issue 
is the determination of the near-crack tip fields for single crystals with several 
active slip systems. These issues, as well as the question of the crack growth, 
will be addressed in our forthcoming papers. Last but not least, it would 
be quite convincing if this theoretical result for the dislocation distribution 
could be compared with the experimental observations and measurements. 
It is hoped that this paper would serve as motivation for experimentalists 
using EBSD-technique to measure the dislocation density near the crack tip. 
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